SPECIFICATION 



Electronic Version 1 .2.8 
Stylesheet Version 1 .0 

Finite-Difference Solver Based on 
Field Programmable Interconnect 

Devices 

Federal Research Statement 

The conditions under which this invention was made are such as to entitle the 
Government of the United States under paragraph 1(a) of Executive Order 1 0096, as 
represented by the Secretary of the Air Force, to the entire right, title and interest 
therein, including foreign rights. 

Background of Invention 

[0001] The present invention is in the field of analog computation circuits, and in 
particular relates to the use of the parasitic resistance of field programmable 
interconnect devices to solve finite difference method problems. 

[0002] The field programmable interconnect device (FPID) is a special-purpose integrated 
circuit, consisting of a large number of transistor-based electronic switches. The FPID 
is generically shown in FIG. 1 . Its design normally involves a number of externally 
available input/output (I/O) terminal contacts, a set of wiring pathways, switches 
between the pathways (represented as circles at a number of the crossing points), and 
a control circuit that determines which switches are closed based a prescribed pattern, 
specified from a configuration port. The FPID permits the flexible and agile 
interconnection between a number of the device's input/output terminals, so that 
normally isolated parts of a networks can be shorted together, or conversely, so that 
in designs, some of the connected parts of a network can be isolated by opening 
switches. 

[0003] 

In the simplified FIG. 1 representation, each terminal is connected to a row and 
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column in the wiring array, so that "A" is actually connected to row A and column A, 
etc. For /7 terminals, this results in 2 //wires (n rows and n columns). Though this 

2 2 
arrangement results in n junctions of rows and columns, only ( n /2- n) switches 

are required to fully connect the n terminals in any combination. This configuration is 
sometimes referred to as a fully connected crossbar. 

[0004] In the unachievable ideal case, the switches represent zero-ohm, zero-length 

wires when closed and infinite resistance connections when opened. Since most FPIDs 
are based on silicon MOSFET devices, however, the switches do not achieve the ideal 
behavior. FIG. 2 illustrates the various representations of the switch in an FPID. FIG. 2a 
is the simplified symbolic representation. FIG. 2b is the familiar standard symbolic 
representation. Figures 2c and 2d represent the n-channel MOSFET and CMOS (n- 
channel plus p-channel) transmission gate structures respectively, which closely 
represent the actual switch structures in FPID devices, FIG. 3 provides a more 
physically accurate representation of an n-channel MOSFET. FIG. 3b illustrates the 
formation of an inversion channel between the drain and source, resulting in a 
conductive path, the situation more closely representing the closure of an FPID switch. 
The switch actually behaves more like a resistor, a fact very important to the principle 
behind the present invention. 

[0005] Since the switch is a poor switch, the FPID is considered a digital device, for use in 
switched logic systems. Switch logic systems compensate for the slight signal 
degradations of transmission gates due to the restorative nature of digital logic 
systems such as CMOS. For general purpose analog, however, the non-zero resistance 
of the transmission gate switch (values may range from 50 ohms to 500 ohms, based 
on the underlying switch design and process technology) results in unwanted signal 
deterioration and design complexity. Hence, even though it is possible to use FPIDs 
for analog applications, it is uncommon to use them for these applications due to the 
normally undesired parasitic resistance. 

[0006] 

It is conceivable, however, that the parasitic resistance could be harnessed in 
particular circuit designs. One such possibility includes the utilization of FPIDs to form 
certain types of resistive networks, in which the normally parasitic resistance now 
plays a key role in the operation of that network. One such circuit class is a linear 
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equation solver, for example, based on the finite difference method. 

[0007] The finite difference method uses a discrete approximation of differential 
equations to reduce them to a system of algebraic equations. For example, the 
following is a derivation of the finite difference representations of Laplace's equation 
in one-dimension. 

[0008] Define Laplace's equation: 

V 2 7 = 0 (i) 

[0009] 
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In one-dimension, Equation (1) becomes: 

The finite forward difference is an approximation of the definition of a derivative: 

dVjx) _ V(x+A)-V(x) = V(x + A)-V(x) 
dx " (x+A)~x A 

Also, define: 

dV(x-A) „ V(x)-V(x-A) (4) 
tic A 

Finite difference representation of higher-level derivatives can be analogously defined: 

dV(x + A) dV(x) 
d 2 V (x) „ dx dx (5) 

dx 2 x 

As A -> 0 ( the approximation improves, being identical to the "true" derivative in the limit 
as A-»0 f Hence, 

dfr dx dx* dx 

etc. So, for convenience we develop the finite difference representation of 

d 2 V(x-A) 
dx 2 

and recognize it as an approximation of 

d 2 V{x) 
dx 2 

using (3) and (4): 

dV{x) dV(x - A) V(x + A)-V(x) V(x)-V(x-A) 

d 2 V(x-A) „ dx dx = A A 

<fr 2 x A 



[0010] 
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V(x + A) + V(x - A) - 2V(x) d 2 V(x - A) d 2 V(x) 

A 2 ~ dx 2 ~ dx 2 {) 

and therefore we can write a finite difference approximation to (1) ( using (6) ) as: 

VT=~ - — ' ^- = 0=>F(x + A) + F(x-A)-2F(x) = 0 

=> F(x + A) + F(x - A) = 2F(jc) 
K(i) B i(F(x + A)-F(x-A)) (7) 

Equation (7) is then the finite difference representation of Laplace's equation in one- 
dimension. 

To simplify implementation in a discrete system or a computer, the A's are typically 
replaced by integral indices, yielding the familiar form of a finite difference equation: 

V(x) = j(V(x+i)-V(x-i)) (8) 

Extending this analysis to multiple dimensions is straightforward. For two dimensions, 
Equation (1) becomes: 

^)-%^^ (9) 

ox oy 

Through the previous analyses, we can directly write the approximation of Equation (9) 
as: 



( V{x + &,y)-V{x,y ) V(x,y)-V{x-^y) \ { ( V{x,y + A)-F(^) V(x,y)-V{x,y-1$ 

V>V(x,y)A * * 11 * * 

A 

_ V(x + &,y) + V(x-A,y) + V(x,y + &) + V{x,y-A)-4y(x,y) 
A 1 

This results in the two-dimensional finite difference method expression: 



V{x,y) = \(V(x + A,y) + V{x - A,y) + V(x,y+A)+ V(x,y- A)) (10) 
4 



[0011] 
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Which can be written in the indexed form as: 



V(x iy ) = -(V(x + l,y) + V(x~hy) + V{x,y + \) + V(x >y -ty (11) 

In particular, it will be shown that it is possible to reduce the solution of Poisson's 
equation: 

V 2 V = p 

over a two-dimensional (2-D) space to an equation at each of many discrete points on a 
grid formed onto this space: 

V (x,y) = 1/4* ( V(x+dx,y) + V(x-dx,y) + V(x,y+dy) + V(x, y-dy) ) 

An electrical analog of this discretization can then be realized by using a grid network of 
resistors. Boundary conditions are simulated by impressing voltages on particular 
nodes. These sources correspond to Dirichlet boundary conditions. Of course, 
Poisson's equation reduces to Laplace's equation in source-free regions: 

V 2 F = 0 

Summary of Invention 

[001 2] In a preferred embodiment, the invention exploits the parasitic resistances of field 
programmable interconnect devices in the form of a programmable resistive grid to 
solve a wide variety of linear partial differential equations. The grid can be 
programmed to mimic the nodal relationships defined in finite difference method 
models with voltages impressed on externally accessible pins corresponding to 
Dirichlet boundary conditions and a means to read out the solutions (voltages) at the 
grid nodal points. A resistive grid may contain up to hundreds of terminals. Problems 
requiring even greater nodal points can be solved sequentially using a plurality of 
resistive grids with the outputs of the first resistive grid component forming the input, 
boundary conditions of the next resistive grid component. Such an approach has a 
distinct advantage over custom solvers that are normally higher in performance but 
fixed in their connection topology. 

[0013] 

Other aspects and advantages of the present invention will become apparent from 
the following detailed description, taken in conjunction with the accompanying 
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drawing, illustrating by way of example the principles of the invention. 

Brief Description of Drawings 

[0014] FIG. 1 illustrates a generic field programmable interconnect device (FPID). 

[001 5] FIG. 2 illustrates various representations of a switch in a FPID, with 2a being a 
simplified symbolic representation, 2b is the standard symbolic representation, 2c 
represents the n-channel MOSFET transmission gate structure, and 2d is the n- 
channel plus p-channel CMOS transmission gate structure. 

[001 6] FIG. 3a illustrates an n-channel MOSFET; 3b shows the formation of an inversion 
channel between the drain (D) and source (S), and 3c illustrates the parasitic 
resistance of the structure. 

[001 7] FIG. 4 is a linear array of connected resistors that is shown to approximate the 
finite difference representation of Laplace's equation in one dimension. 

[001 8] FIG. 5a is a cross-sectional representation of a capacitor. 

[001 9] FIG. 5b is a model of the FIG. 5a capacitor as a resistive network. 

[0020] FIG. 6a is a two-dimensional mesh resistor that approximates the finite difference 
representation of Laplace's equation in two dimensions. 

[0021] FIG. 6b is a north-east-west-south (NEWS) resistor grid pattern of nearest 
neighbors. 

[0022] FIG. 7a is an electrical analog using a grid network of resistors to represent the 
solution of Poisson's equation over a 2-D space reduced to discrete points on a grid. 

[0023] FIG. 7b depicts a smaller sub-region of FIG. 7a. 

[0024] FIG. 7c illustrates notionally how switches might be set in a FPID switch array to 
produce the FIG, 7b sub-region. 

[0025] FIG. 8a illustrates a small sample network in which a number of nodes are forced 
with a V=0 Dirichlet boundary condition and a single node is forced with a V=1 
Dirichlet boundary condition. 
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[0026] FIG. 8b is a table summarizing the required switch closures for the FPID network 
shown in FIG. 8a. 

[0027] FIG. 9a represents a 2-D confined charge box sample problem with the box 
shown in cross-section. 

[0028] FIG. 9b is a table of the pin connections and expected measurements for the 9a 
problem solved by using a FPID device. 

[0029] FIG. 1 0 illustrates a self-contained FPID with a number of digital-to-analog 

converters and analog-to-digital converters added to the FPID to produce a hybrid 
computer. 

[0030] FIG. 1 1 illustrates temporary connections, made one-at-a-time, to acquire voltage 
measurements for a FPID. 

[0031] FIG. 12 illustrates a larger network comprised of a plurality of individual FPID's 
interconnected with reconfigurable terminal pins. 

[0032] FIG. 1 3 shows a finite difference mesh (1 1 a) that is too large to fit onto one FPID 
device and 1 1 b shows this mesh might be sectioned into four segments. 

[0033] FIG. 1 4 shows the solution process for the first complete cycle of a four segment 
finite difference mesh, where 1 4a is the first segment corresponding to FIG. 1 3a, 1 4b 
corresponds to 1 3b, etc. 

Detailed Description 

The purpose of the current invention is to exploit the undesired parasitic 
resistance in a special class of digital integrated circuits to perform analog 
computation. This class of integrated circuits, referred to as field programmable 
interconnect devices (FPID's) act as crossbar switches, permitting the arbitrary 
connection of many signals attached to the package pins of the components. The 
invention extends the application of this crossbar to computation, which can now be 
used as a building block in analog or hybrid (analog plus digital) computer 
architectures. 

The present invention exploits the parasitic resistance of FPIDs to form a 
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programmable resistive grid that can be used to solve certain finite difference method 
(FDM) problems. A FPID component contains dozens to hundreds of terminals that can 
be shorted together or isolated under program control. "Shorting" and "isolating" are 
relative concepts. Since contemporary FPID components are built using MOSFET 
transistors in silicon integrated circuits configured as switches, their connective paths, 
formed by path closures, have intrinsic channel resistance. So, instead of "shorting" 
terminals together, the path closure forms a resistance. This resistance is generally 
undesirable, but as it is nominally consistent from switch to switch, it is possible to 
exploit this regularity to form resistive grids. The grids take on configurations as 
defined by programming the connection paths between any given combination of 
terminals. Similarly, when the path is opened, an extremely large but finite resistance 
remains in effect between terminals. Fortunately, for the purposes of the current 
invention, it is possible to neglect this effect, i.e., to treat the open condition as the 
ideal case of infinite resistance or isolation. 

[0036] Equation 7 was previously shown to be the finite difference representation of 
Laplace's equation in one dimension. It will now be shown that a linear array of 
connected resistors implements the approximation of Laplace's equation given in Eq. 
7.The network of resistors is defined in FIG. 4. Next, we write equations for V(x) and V 
(x+ A ) based on elementary circuit theory: 

V(x) = V(x-A) + !R 

Combining these equations yields Equation (7): 

V(x + A) = V(x) + (V(x) - V(x - A)) F(jc) = ^(V(x + A) + V{x - A)) 

The A's are typically replaced by integral indices, yielding the familiar form of a finite 
difference equation: 

V(x)A(V(x+i)-V(x-i)) (8) 

[0037] Consider a very simple 1 -D exemplary problem, i.e., solving the electrostatic 

potential in an infinite slab. In this case, two conducting slabs are provided to contain 
a dielectric slab, one positioned at x=0 and one at x=l . Though shown as finite, the 
slabs are understood to have infinite extent in the y and z axes (an ideal parallel plate 
capacitor). Boundary conditions are provided in the form of specified voltages at each 
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conducting slab, namely V(0)=0 and V(l)=l . This is shown in FIG 5a. The 
corresponding differential equation problem is simply: 



with V(0) = 0 and V(l) = \. 

An analytic solution is quickly developed by integrating this equation, producing 

and once again integrating to produce 

V(x)= ^k l dx = k l x + k 2 

it is obvious that V(0)=0 implies that k 2 =0, and V(1)=1 implies that ki=1, producing the 
very simple result that 

V(x) = x becomes a reference verification for a finite difference model representation of 
the initial problem. Based on the previous discussion, it is straightforward to model the 
slab as a resistive network. 

In this case, the resistances represent discretization of the slab at nine 
equidistant points as shown in FIG. 5b. The voltage applied to point 1 equals zero volts, 
whereas the voltage applied to point nine equals one volt. Whether a discrete finite 
difference model is used in which 

[0038] or the resistor network is built and measured or analyzed, clearly the resulting 
findings will be as follows from points 1 to 9: 0, 0.125, 0.250, 0.375, 0.500, 0.625, 
0.750, 0.875, 1 ,000. This corresponds exactly to the analytic solution previously 
specified by V(x) = x. 

[0039] For the two-dimensional (2-D) case, a 2-D mesh resistor (FIG. 6a) implements the 
finite difference representation for Laplace's equation in 2-D, previously derived in 
index form as: 

K(^^)=I(K(x+l ? >;) + K(x-l ) j) + F(x^ + l) + F(x s y-l)) (11) 
4 

[0040] 

As such, a FDM grid can be imposed by establishing a north-east-west-south 
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(NEWS) pattern of nearest neighbors starting with a chosen terminal (see FIG. 6b) until 
all the terminals are consumed by the network formed. This resistive network models 
a Laplacian 2-D partial differential equation. More elaborate, multi-dimensional (when 
viewed from a Euclidean perspective) grids can easily be established by simply 
defining the appropriate connective paths. 

[0041] The present invention then deliberately exploits the parasitic resistance of FPID's 
to form a partial differential equation solver. This involves configuring the FPID to 
create networks similar to that shown in FIG. 7 by closing the appropriate set of 
switches. By closing FPID switches in this manner, a network of parasitic resistances 
results. While the parasitic resistance values may vary from one design to another, or 
even from one lot run to another, the parasitic resistance values across a particular 
device typically are very closely matched due to process controls used to maximize 
device yield. Therefore, the FPID can be used as an effective analog-domain equation 
solver, which can converge much faster than a corresponding digital implementation 
of the finite difference method. 

[0042] FIG. 7a illustrates the resistive grid topology typical of 2-D discretized 

implementations of a finite difference method equation that can be implemented 
using the parasitic resistances of switches in a FPID. FIG. 7b depicts a smaller sub- 
region, and FIG. 7c illustrates notionaliy how switches might be set in a FPID switch 
array to produce the FIG. 7b sub-region. 

[0043] The implementation of the FPID equation solver is illustrated simply with a small 
example. FIG. 8a shows a network in which a number of nodes are forced with a V=0 
Dirichlet boundary condition and a single node is forced with a V=1 Dirichlet 
boundary condition. Two exterior pins are required, a ground pin connected to 
terminal number 1 7 for the V=0 boundary condition and a connection of a one volt 
source to terminal number 4 for the V=l boundary condition. No other external 
connections are required to the FPID, except for the connections required to power 
the device control circuitry and the configuration port. The pin connections between 
nodes are then prescribed in an obvious way. For example, pin number 6 is connected 
to pins number 2, 5, 7, and 10. Any ground connection is achieved by connecting to 
pin number 1 7. The table (FIG. 8b) summarizes the required switch closures. 
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[0044] Once these switch closures have been accomplished, the FIG. 8 network is formed 
and the solution of Laplace's / Poisson's equation is available very quickly (much less 
than one microsecond) and can be directly evaluated at any terminal through 
measurement of the voltage at particular pins. Changing the boundary conditions can 
be done continuously and independently, without changing the switch configuration. 
It is also possible to completely alter the network connections, representing a 
different equation and boundary conditions. 

[0045] While the 1 -D capacitor case above demonstrates a nearly trivial example, it is 
equally straightforward to extend the same technique to much more analytically 
difficult examples. FIG. 9a shows a relatively simple two-dimensional confined charge 
box. A cross-sectional view of an infinitely long box is shown. An analytic solution of 
V(x,y) within the box is not easily developed. However, a discretized version can be 
solved by conventional techniques using a matrix of linear equations or a discrete 
finite difference numerical method. The problem can also be much more rapidly 
solved directly using a FPID device. The pin connection and expected measurements 
are provided in the FIG. 9b table. 

[0046] In this case, the pin number corresponds not necessarily to the identically 

numbered pin on the FPID device, but rather to a specific sequence on sixteen user 
pins. The voltages on pin #1 and pin #7 (FIG. 9a) are fixed as boundary conditions, 
and the voltages on other pins are measured after the FPID is suitably programmed to 
model the neighborhood relationships implied in the shown discretization. 

[0047] It is important to indicate that there are practical limits to the use of the FPID as 
an analog-domain equation solver. First, the excursion range or operational window 
of an FPID is limited. Therefore the range of voltages impressed upon the FPID must 
not exceed in a positive or negative polarity the magnitude that would cause the 
transistor switches in an FPID to breakdown the gate voltage or forward bias the 
substrate connection, for example. A typical industrial FPID device might permit a 
voltage range from 0V to IV for example. 

[0048] j 0 S j m p|jfy anc j automate the creation and instrumentation of test partial 

differential equation configurations, it is possible to make a self-contained circuit, 
such as shown in FIG. 10. Here, a number of digital-to-analog converters (DACs) and 
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analog-to-digital converters (ADCs) are added to the FPID to produce a hybrid 
computer. A number of obvious design details are omitted, such as: (1) the bussing 
arrangements and enable signals for the different peripheral blocks, (2) that the 
operating windows of the ADC/DAC components must be aligned with the effective 
operating voltage range of the FPID switch, and (3) the DAC outputs must be 
themselves isolated (using possible a transmission gate switch) when they are not 
used actively in a particular problem formulation. The primary function of the DACs is 
to inject boundary conditions, while the function of the ADCs is to read analog 
voltages corresponding to nodal solutions. Therefore, the number of nodal boundary 
conditions is limited by the number of DAC circuits. The number of ADCs may be as 
few as one. If the ADC employs an infinite impedance front end, then it is conceivable 
that an ADC can be tied to a single pin of the FPID. It is then possible, one by one, to 
close switches between the measurement node (defined as the single pin connected to 
an analog-to-digital converter) and particular nodes in an active equation under 
solution for the purposes of measurement. This configuration is suggested in FIG. 1 1 
where temporary connections are made one at a time to acquire voltage 
measurements. When the input impedance 18 of the ADC is infinite, then the parasitic 
resistance formed by the temporarily switched measurement connections will not 
contribute to measurement error. In the case of finite resistance, a systematic 
(correctable) error due to the corresponding voltage ladder effects. 

[0049] It is sufficient, therefore, to obtain all measurements with one ADC per FPID. The 
addition of more ADCs is only advantageous in the cases where it is necessary to 
more rapidly acquire signal measurements. 

[0050] | t js p 0SS jb| e t; 0 extend these concepts by adding FPIDs. In this case, the external 
(reconfigurable) terminal pins of FPIDs are interconnected in some way to facilitate the 
extension of the parasitic networks of the ensemble to a larger effective network. This 
possibility is suggested in FIG. 1 2. Two notes are pointed out here. First, the 
configuration by which the FPIDs are interconnected is shown as a 2-D planar mesh. 
In fact, an almost arbitrary number of arrangements are possible, including those 
configurations where pins are shared with more than two FPIDs. Second, the 
configuration portals of each of the FPIDs are shown as independent. They must be 
connected to a computer or programming source to supply the commands for switch 
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closures within the FPIDs. It is possible that, rather than distinct and independent, that 
the configuration ports might be bussed together or daisy-chain connected together, 
consistent with the practices used in complex systems containing multiple FPCA 
devices. 

[0051] The primary motivation for using multiple FPIDs is to extend the size of the 

problem that could be solved using the analog-domain approach that is the basis of 
the invention. It is also possible to extend the size of solvable problems by extending 
the FIG. 7 configuration, using a DAC to optionally drive every terminal pin of an FPID. 
To illustrate how such an approach can be used to solve an "oversized" problem, a 
simple example is provided in FIG. 13. In this case, a finite difference mesh (FIG, 13a) 
is too large to fit onto one FPID device. Instead, the problem is sectioned into four 
pieces or "segments", as suggested in FIG. 13b. The first steps of an iterative solution 
process are shown in FIG. 14. Briefly, the solution process involves driving the 
boundaries of the current "segment" with values obtained from simulating the 
adjacent segments in the last cycle. The DACs are used to force these values 
temporarily during a given cycle. These values will typically change each cycle, and so 
they are called floating (Dirichlet) boundary conditions. Initial conditions may be 
chosen empirically or simply set to null (undriven). 

[0052] In FIG. 14, the first complete cycle of a four-segment simulation is shown. The 
first segment, corresponding to the upper left quadrant of FIG, 1 3b, is configured or 
loaded onto the FPID, including any external boundary conditions, which corresponds 
to FIG. 14a. The FPID is operated very briefly and the boundary cell values are 
measured and stored. The next piece, corresponding to the upper right quadrant of 
FIG. 1 3b, is loaded onto the FPID as suggested in FIG. 14b. In this case, some values 
from the FIG. 14a simulation, which were measured, are now used themselves as 
boundary conditions to be driven by the appropriate DACs, providing a more 
approximate boundary condition set for the new simulation. This new simulation in 
turn produces node values to be measured and used in subsequent simulations. The 
process is repeated in FIGURES 14c and 14d, corresponding to the lower left and right 
quadrants, respectively, of FIG. 1 3b. In subsequent cycles of the computation process 
illustrated in FIG, 14, the boundary measurements eventually converge to accurate 
values for the overall simulation. 
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